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We obtain the finite-temperature unconditional master equation of the density matrix for two 
coupled quantum dots (CQD) when one dot is subjected to a measurement of its electron occupation 
number using a point contact (PC) . To determine how the CQD system state depends on the actual 
current through the PC device, we use the so-called quantum trajectory method to derive the zero- 
temperature conditional master equation. We first treat the electron tunneling through the PC 

in : 

barrier as a classical stochastic point process (a quantum-jump model). Then we show explicitly 
that our results can be extended to the quantum-diffusive limit when the average electron tunneling 
rate is very large compared to the extra change of the tunneling rate due to the presence of the 
electron in the dot closer to the PC. We find that in both quantum-jump and quantum-diffusive 
cases, the conditional dynamics of the CQD system can be described by the stochastic Schrodinger 
, equations for its conditioned state vector if and only if the information carried away from the CQD 

^ ■ system by the PC reservoirs can be recovered by the perfect detection of the measurements. 

a 

a 

I. INTRODUCTION 

The origins and mechanisms of decoherence (dephasing) for quantum systems in condensed matter n nhysics have 
attracted much attention recently due, to a number of studies in nanostructure mesoscopic systemsEHj and vari- 
ous proposals-fpr. quantum computersoH. One of the issues is the connection., between decoherence and quantum 
measurementllj'Ell for a quantum system. It was reported in a recent experiment with a "which-path" interferometer 
that Aharonov-Bohm interference is suppressed owing to the measurement of which path an electron takes through 
the double-path interferometer. A biased quantum point contact (QPC) located close to a quantum dot, which is 
built in one of the interferometer's arms, acts as a measurement device. The change of transmission coefficient of the 
■ QPC, which depends on the electron charge state of the quantum dot, can bepdetepted. The decoherence rate due to 
\ the measurement by the QPC in this experiment has been calculated in Refs.ll3il3. 

A quantum-mechanical two-state system, coupled to a dissipative environment, provides a universal model for many 
physical systems. The indication of quantum coherence can be regarded as the oscillation or the interference between 
the probability amplitudes of finding a particle between the two states. In this paper, we consider the problem of 
an electron tunneling between two coupled quantum dots (CQDs) using a low-transparency point contact (PC) or 
tunnel junction as a detector (^environment) measuring the position of the electron (see Fig. Q). This problem has been 
Q extensively studied in Refs.li3 ciL |X he . ca se of measurements by a general QPC detector with arbitrary transparency has 
also been investigated in RefsJl3il3rprj?l .^pad^tiotL-a similar system measured by a single electron transistor rather 
than a PC has been studied in RefsEOOE^Gj'EffEaEj. The influence of the detector (environment) on the measured 
system can be determined by the reduced density matrix obtained by tracing out the environmental degrees of the 
freedom in the total, system plus environment, dejasity matrix. The master equation (or rate equations) for this CQD 
system have been derived and analyzed in RefsJIS'ta (here we refer to the rate equations as the first order differential 
equations in time for both diagonal and off-diagonal reduced density matrix elements). This (unconditional) master 
equation is obtained when the results of all measurement records (electron current records in this case) are completely 
ignored or averaged over, and describes only the ensemble average property for the CQD system. However, if a 
measurement is made on the system and the results are available, the state or density matrix is a conditional state 
conditioned on the measurement results. Hence the deterministic, unconditional master equation cannot describe 
the conditional dynamics of the CQD system in a single realization of continuous measurements which reflects the 
stochastic nature of an electron tunneling through the PC barrier. Consequently, the conditional master equation 
should be employed. In condensed matter physics usually many identical quantum systems are prepared at the 
same time and a measurement is made upon the systems. For example, in nuclear or electron magnetic resonance 
experiments, generally an ensemble of systems of nuclei and electrons are probed to obtain the resonance signals. This 
implies that the measurement result in this case is an average response of the ensemble systems. On the other hand, 
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for various proposed condensed-matter quantum computer architecturesETtl, how to readout physical properties of a 
single electronic qubit, such as charge or spin at a single electron level, is demanding. This is a non-trivial problem 
since it involves an individual quantum particle measured by a practical detector in a realistic environment. It is 
particularly important to take account of the decoherence introduced by the measurements on the qubit as well as 
to understand how the quantum state of the qubit, conditioned on a particular single realization of measurement, 
evolves in ti me for the purpose of quantum computing. 

KorotkovtSEil has obtained the Langevin rate equations for the CQD system. These rate equations describe the 
random evolution of the density matrix that both conditions, and is conditioned by, the PC detector output. In 
his approach, the individual electrons tunneling through the PC barrier were ignored and the tunneling current was 
treated as a continuous, diffusive variable. More precisely, he considered the change of the output current average 
over some small time r, (I), with respect to the average current Ii, as a Gaussian white noise distribution. He then 
updated (J) in the density-matrix elements using the new values of (/) after each time interval r. However, treating 
the tunneling current as a continuous, diffusive variable is valid only when the average electron tunneling rate is very 
large compared to the extra change of the tunneling rate due to the presence of the electron in the dot closer to 
the PC. The resulting derivation of the stochastic rate equations is semi-phenomenological, based on basic physical 
reasoning to deduce the properties of the density matrix elements, rather than microscopic. __ . .. . 

To make contact with the measurement output, in this paper we present a quantun^KaiectoryE^^TCorEB measurement 
analysis to the CQD system. We first use the quantum open system approachEirEil to obtain the unconditional 
Markovian master equation for the CQD system, taking into account the finite-temperature effect of the PC reservoirs. 
Particularly, we assume the transparency of the PC detector is small, in the tunnel-junction limit. Subsequently, we 
derive microscopically the zero-temperature conditional master equation by treating the electron tunneling through 
the PC as a classical stochastic point process (also called a quantum-jump model)EZrc3€3. Generally the evokttion of 
the system state undergoing quantum jumps (or other stochastic processes) is known as a quantum trajectory^. Real 
measurements (for example the photon number detection) that correspond approximately to the ideal quantum-jump 
(or point-process) measurement are made regularly in experimental quantum optics. For almost all-infinitesimal time 
intervals, the measurement result is null (no photon detected). The system in this case changes infinitesimally, but 
not unitarily. The nonunitary component reflects the changing probabilities for future events conditioned on past 
null events. At randomly determined times (conditionally Poisson distributed), there is a detection result. When this 
occurs, the system undergoes a finite evolution, called a quantum jump. In reality these point processes are not seen 
exactly due to a finite frequency response of the circuit that averages each event over some time. Nevertheless, we 
first take the zero-response time limit and consider the electron tunneling current consisting of a sequence of random 
S function pulses, i.e., a series of stochastic point processes. Then we show explicitly that |0«£-|results can be extended 
to the quantum-diffusive limit and reproduce the rate equations obtained by KorotkovlialHd. We refer to the case 
studied by Korotkovllac3 as quantum diffusion, in contrast to the case ofrquantum jumps considered her e ,.. H ence 
our quantum trajectory approach may be considered as a formal derivations of the rate equations in Refs.EaE3. Wc 
find in both quantum-jump and quantum-diffusive cases pti^|-the||Gonditional dynamics of the CQD system can be 
described by the stochastic Schrodinger equations (SSEs)OEJ'EZIoc3 for the conditioned state vector, provided that 
the information carried away from the CQD system by the PC reservoirs can be recovered by the perfect detection of 
the measurements. 

This paper is organized as follows. In Sec. 0, we sketch the derivation of the finite-temperature unconditional master 
equation for the QCD system. To determine how the CQD system state depends on the actual current through the 
PC device, we derive in Sec. Ill the zero-temperature conditional master equation and the SSE in the quantum- 
jump model. Then in Sec. [iv| we extend the results to the case of quantum diffusion and obtain the corresponding 
conditional master equation and SSE. The analytical results in terms of Bloch sphere variables for the conditdjmal 
dynamics are presented in Sec. [v| Specifically, we analyze in this section the localization rate and mixing ratezj'EiJcj. 
Finally, a short conclusion is given in Sec. VI. Appendix His devoted to the demonstration of the equivalence between 

rka and t 



the conditional stochastic rate equations in Refs.IL2rei!l and those derived microscopically in the present paper. 



II. UNCONDITIONAL MASTER EQUATION FOR THE CQD AND PC MODEL 



The appropriate way to approach quantum measurement problems is to treat the m e asure d system, the detector 
(environment), and the coupling between them microscopically. Following from Refs.c3E3iH3, we describe the whole 
system (see Fig. [l]) by the following Hamiltonian: 

H = Ti-CQD + Ti-PC + T~tcoup (1) 

where 
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'HCQD = ft [wicjci + 0>2<4 C 2 + ^(cl c 2 + 4 c l)] i ( 2 ) 
Hpc = ft a Lk a Lk + ^k a Rk a Rk) + ( T kq a Lk a nq + T qk a R q a Lk) , (3) 



k k.q 

H coup = c tci (xk q a[ k a Rq + X* qk a Rq a Lk) ■ (4) 

k,q 

7~Lcqd represents the effective tunneling Hamiltonian for the measured CQD system. For simplicity, we assume strong 
inner and inter dot Coulomb repulsion, so only one electron can occupy this CQD system. We label each dot with an 
index 1,2 (see Fig. |l|) and let Cj (cj) and Tvjii represent the electron annihilation (creation) operator and energy for 
a single electron state in each dot respectively. The coupling between these two dots is given by Ml. The tunneling 
Hamiltonian for the PC detector is represented by TLpc where a^k-, dRk an d %UJ^, Uoj^ are respectively the fermion 
(electron) field annihilation operators and energies for the left and right reservoir states at wave number k. One 
should not be confused by the electron in the CQD with the electrons in the PC reservoirs. The tunneling matrix 
element between states k and q in left and right reservoir respectively is given by T kq . Eq. (f|), 7i COMp , describes the 
interaction between the detector and the measured system, depending on which dot is occupied. When the electron 
in the CQD system is close by to the PC (i.e., dot 1 is occupied), there is a change in the PC tunneling barrier. This 
barrier change results in a change of the effective tunneling amplitude from T kq — > T kq + Xkq- As a consequence, 
the current through the PC is also modified. This changed current can be detected, and thus a measurement of the 
location of the electron in the CQD system is effected. 

The total density operator R(t) for the entire system in the interaction picture satisfies: 

i _ 1 



Ri(t) = — [Hi(t), Ri(0)} - ^ I dt' [H!(t), [Hj(t , ),R I (t / )}} . (5) 



The dynamics of the entire system is determined by the time-dependent HamiltonianB 



H x (t) = ( T ki + XkAa) aifcW^-O* + H.C., (6) 

k,q 

where we have treated the sum of the tunneling Hamiltonian parts in TLpc and Ti. coup as the interaction Hamiltonian 
Ttj, and H.C. stands for Hermitian conjugate of the entire previous term. By tracing both sides of Eq. (|^) ovepthe. 
bath (reservoir) variables and then changing from the interacting picture to the Schrodinger picture, we obtainEllLJ 
the finite-temperature, Markovian master equation for the CQD system: 

p{t) = -^{H C QD,p(t)} + V[T + + X +ni ]p{t) + V[T* + Af*m]/j(t), (7) 

where p(t) = Trp,R(t) and Tr^ indicates a trace over reservoir variables. In arriving at Eq. (Q), we have made the 
following assumptions and approximations: (a)treating the left and right fermion reservoirs in the PC as thermal 
equilibrium free electron baths, (b)weak system-bath coupling, (c)small transparency of the PC, i.e., in the tunnel- 
junction limit, (d)uncorrclatcd and factorizable system-bath initial condition (e)relaxation time scales of the reservoirs 
being much shorter than that of the system state, (f)Markovian approximation, (g)|eV|, kpT <C Pl(r)i and (^energy- 
independent electron tunneling amplitudes and density of states over the bandwidth of max(|eV^|, ksT). Here ks is 
the Boltzmann constant, T represents the temperature, eV = fit — pr is the external bias applied across the PC, and 
Pl and pr stand for the chemical potentials in the left and right reservoirs respectively. In Eq. (Q), n\ = c\c\ is the 
occupation number operator for dot 1. The parameters T± and X± are given by 

\T±\ 2 = D ± = 2ne\T m \ 2 g L g R V ± /h, (8a) 
\T± +X ± \ 2 = D' ± = 27re|T o + Xoo\ 2 gLg R V ± /h, (8b) 

where D± and D' ± are the average electron tunneling rates through the PC barrier in positive and negative bias 
directions at finite temperatures, without and with the presence of the electron in dot 1 respectively. Here the 
effective finite-temperature external bias potential, eV± is given by the following expression: 

CV± = l-exp[ T eV/(k B T)Y (9) 

Too and xoo are energy-independent tunneling amplitudes near the average chemical potential, and and gR are 
the energy-independent density of states for the left and right fermion baths. Note that the average electron currents 
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through the PC barrier is proportional to the difference between the average electron tunneling rate in opposite 
directions. Hence, the avawtge currents eD = e(D + — D ) and eD' = e(D' + — D'_), following from Eq. (|J) and (^), are 
temperature independentc3c3 at least for a range of low temperatures ksT <C Pl{r)- ln addition, the current-voltage 
characteristic in the linear response region \eV\ <C Pl(r) is of the same form as for an Ohmic resistor, though the 
nature of charge transport is quite different in both cases. ___ _ 

We have also introduced, in Eq. (Q), an elegant superoperatorElKcZrca T>, widely used in measurement theory in 
quantum optics. Physically the "irreversible" part caused by the influence of the environment in the unconditional 
master equation, is represented by the T> superoperator. Generally superoperators transform one operator into another 
operator. Mathematically, the expression T>[B]p means that superoperator V takes its operator argument B, acting 
on p. Its precise definition is in terms of another two superoperators J and A: 

V[B]p = J[B]p - A[B}p, (10) 

where 

J[B]p = BpB\ (11) 

A[B]p= (B^Bp + P B^B)/2. (12) 

The form of the master equation (Q), defined through the superoperator 2? [£?]/?(£), .-preserves the positivity of the 
density matrix operator p(t). Such a Markovian master equation is called a LindbladEl form. _ 

To demonstrate the equivalence between the master equation^-/|7|) and the rate equations derived in Ref.Ej, we 
evaluate the density matrix operator in the same basis as in Ref.EJj and obtain 

Paa(t) = ifl[pab(t) - Pba(t)} , (13a) 

Pab {t) = ie Pab {t) + m[ Paa {t) - p bb (t)} - {\X T \ 2 /2)p ab {t) + ilm{T*X+ - T*X_) Pab (t) (13b) 

Here he — H(lu2 — oj\) is the energy mismatch between the two dots, Pij{t) = (i\p(t)\j), and p aa (t) and p bb (t) are 
the probabilities of finding the electron in dot 1 and dot 2 respectively. The rate equations for the other two density 
matrix elements can be easily obtained from the relations: p bb (t) = 1 — p aa (t) and p ba (t) — P* ab (t)- Compared to 
an isolated CQD system, the presence of the PC detector introduces two effects to the CQD system. First, the 



imaginary part of (T^X + — T1X ) (the last term in Eq. (13b)) causes an effective temperature-independent shift in 



the energy mismatch between the two dots. Here, (T^_X + — T*X_) — T*X is a temperature-independent quantity, 
where T = T + (0), X = X + (0); i.e., T + and X + evaluated at zero temperature respectively. Second, it generates a 
decoherence (dephasing) rate 

T d = \X T \ 2 /2 (14) 

for the off-diagonal density matrix elements, where \Xt\ 2 = \X + \ 2 + \X-\ 2 . We note that the decoherence rate 
comes entirely from the effect of the measurement revealing where the electron in the CQDs is located. If the 
PC detector does not distinguish which of the dots the electron occupies, i.e., X± = Ojpthen Yd = 0. The rate 
equations in Eq. (|ll|) are exactly the same as the zero-temperature rate equations in Ref£3 if we as sume that the 
tunneling amplitudes are real, T ao — Tq and x o = Xoo- I n that case, the last term in Eq. ( |l3b| ) vanishes and 
= X 2 /2 = [yW — \ /f D) 2 /2. Actually, the relative phase between the two complex tunneling amplitudes may 
produce additional effects on conditional dynamics of the CQD system as well. This will be shown later when we 
discuss conditional dynamics. Physically, the presence of the electron in dot 1 raises the effective tunneling barrier 
of the PC due to electrostatic repulsion. As a consequence, the effective tunneling amplitude becomes lower, i.e., 
D' = \T + X\ 2 < D = \T\ 2 . This sets a condition on the relative phase 9 between X and T : <eos 6 < -\X\/(2\T\). 

The dynamics of the unconditional rate equations at zero temperature was analyzed in Ref.t£l. Here, following from 
Eqs. (^4(), (g) and ([})), we find that the temperature-dependent decoherence rate due to the PC thermal reservoirs 
has the following expression: 



r d (0) eV \2k B T 

As expected, Td{T) increases .with increasing temperature, although the average tunneling current through the PC is 
temperature independenta'C^I for the same range of low temperatures fc^T <C Pl(r)- This temperature dependence 
of the decoherence rate is in fact just the temperature dependence-of the zero-frequency noise power spectrum of 
the current fluctuation in a low-transparency PC or tunnel junctionEil. The CQD system weakly coupled to another 
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finite-temperature environment beside the PC detector was discussed in Ref.E3. However, the influence of the finite- 
temperature PC reservoirs on the CQD system, presented here, was not taken into account. The finite-temperature 



decoherence rate jOi an one^e 
calculated in Ref.EJ. In RefH 



ectron state in a quantum dot due to charge fluctuation of a general QPC has been 
, the temperature-dependent decoherence rate for a two-state system caused by a QPC 



detector has been discussed specifically in the context of the measurement problem. 



III. QUANTUM-JUMP, CONDITIONAL MASTER EQUATION 

So far we have considered the evolution of reduced density matrix when all the measurement results are ignored, 
or averaged over. To make contact with a single realization of the measurement records and study the stochastic 
evolution of the quantum state, conditioned on a particular measurement realization, we derive in this section the 
quantum-jump, conditional master equation at zero temperature. 

The nature of the measurable quantities, such as accumulated number of electrons tunneling through the PC barrier, 
is stochastic. On average of course the same current flows in both reservoirs. However, the current is actually made 
up of contributions from random pulses in each reservoir, which do not necessarily occur at the same time. They 
are indeed separated in time by the times at which the electrons tunnel through the PC. In this section, we treat 
the electron tunneling current consisting of a sequence of random S func tion , pulses. In other words, the measured 
current is regarded as a series of point processes (a quantum-jump model)EHSEl The case of quantum diffusion will 



be analyzed in Sec. [V 



Before going directly to the derivation, we discuss some general ideas concerning quantum measurements. If the 
system under observation is in a pure quantum state at the beginning of the measurement, then it will still be in a 
pure conditional state after the measurement, conditioned on the result, provided no information is lost. For example, 
if the initial normalized state is \ip(t)), the unnormalized final state given the result a at the end of the time interval 
[t, t + dt) of the measurement becomes 

\i> a [t + dt)) = M a {dt)\iP{t)), (16) 
where {M a (t)} represents a set of operators which define the measurements and satisfies the completeness condition 

Y J Ml{t)M a {t) = l. (17) 



Eq. (|17j) is simply a statement of conservation of probability. The corresponding unnormalized density matrix, 
following from Eq. ([l6]), is given by 

p a (t + dt) = \^ a {t + dt))(rp a (t + dt)\ = J[M a (dt)]p(t), (18) 

where p(t) — \tj){t)){^}{t)\ and the superoperator J is defined in Eq. (|TT|). Of course, if the measurement is made 
but the result is ignored, the final state will not be pure but a mixture of the possible outcome weighted by their 
probabilities. Consequently, the unconditional density matrix can be written as 

p(t + dt) = ^p a {t + dt) = ^Pr[a]p Q (< + dt), (19) 

a ol 

where Pr[a] = Tr[p a (t + dt)] stands for the probability for the system to be observed in the state a, and p a (t + dt) = 
p a (t + dt)/Pr[a] is the normalized density matrix at time t + dt. 

Now we proceed to derive the quantum-jump, conditional master equation in the following. Only two measurement 
operators M a (dt) for a — 0, 1 are needed for a measurement record which is a point process. For most of the 
infinitesimal time intervals, the measurement result is a = 0, regarded as a null result. On the other hand, at 
randomly determined times, there is a result a = 1, referred as a detection of an electron tunneling through the PC 
barrier. Formally, we can write the current through the PC as 

i(t)=edN(t)/dt, (20) 

where e is the electronic charge and dN(t) is a classical point process which represents the number (either zero or 
one) of tunneling events seen in an infinitesimal time dt. We can think of dN(t) as the increment in the number of 
electrons N(t) in the drain-in time dt. It is this variable, the accumulated electron number transmitted through the 
PC, which is used in Refs£3'E3E3. The point process is formally defined by the conditions on the classical random 
variable dN c (t): 
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[dN c {t)\ 2 = dN c (t), (21) 
E[dN c (t)] = Tr[p lc (t + dt)] = Tr{J[Mi{dt)]p D (t)} = V lc (t)dt. (22) 

Here we explicitly use the subscript c to indicate that the quantity to which it is attached is conditioned on previous 
measurement results, the occurrences (detection records) of the electrons tunneling through the PC barrier in the 
past. E[Y] denotes an ensemble average of a classical stochastic process Y. Eq. (|2l"l) simply states that dN c (t) equals 
either zero or one, which is why it is called a point process. Eq. ( p^ ) indicates that the ensemble average of dN c (t) 
equals the probability (quantum average) of detecting electrons tunneling through the PC barrier in time dt. In 
addition, dN c (t) is of order dt and obviously all moments (powers) of dN c (t) are of the same ord er as dt. Note here 



that the density matrix p c (t) is not the solution of the unconditional reduced master equation, Eq. ( 25a ). It is actually 
conditioned by dN c (t') for t' < t. 

The stochastic conditional density matrix at a later time t + dt can be written as: 

*« + di) = dNM mu^k + 1 1 - "^row < 23 > 

Eq. (|23|) states that when dN c (t) = (a null result), the system changes infinitesimally via the operator M (dt) and 
hence the p c (t + dt) = po c (t + dt). Conversely, if dN c (t) = 1 (a detection), the system goes through a finite evolution 
induced by the operator M\(dt), called a quantum jump. The corresponding normalized conditional density matrix 
then becomes p\ c {t + dt). One can see, with the help of Eqs. (pp|), that in this approach the instantaneous system 
state conditions the measured current (see Eq. (^||)) while the measured current itself conditions the future evolution 
of the measured system (see Eq. (p3[)) in a self-consistent manner. It is straightforward to show that the ensemble 
average of the conditional density matrix equals the unconditional one, E[p c (t)] — p(t). Tracing over both sides of 
Eq. ( |l9| ) for a = 0, 1, we obtain 

Tr[p 0c (t + dt)} = 1 - Tr[p lc (i + dt)]. (24) 

Then taking the ensemble average over the stochastic variables dN c (t) on both sides of Eq. (^3|), replacing E[dN c (t)] 
by using Eq. ( p2| ) , and comparing the resultant equation with Eq. Jl9| ) completes the verification. 

Next we find the specific expression of pi c {t + dt) and po c (t + dt) and derive the conditional master equation for 
the CQD system measured by the PC. If a perfect PC detector (or efficient measurement) is assumed, then whenever 
an electron tunnels through the barrier, there is a measurement record corresponding to the occurrence of that event; 
there are no 'misses' in the count of the electron number. As a result, the information lost from the system to 
the reservoirs can be recovered using a perfect detector. Here we assume a zero-temperature case for the efficient 
measurement. At finite temperatures, the electrons can, in principle, tunnel through the PC barrier in both directions. 
But experimentally the detector might not be able to detect these electron tunneling processes on both sides of the 
PC barrier. This may result in information loss at finite temperatures. Hence, at zero temperature the unconditional 
master equation (0) reduces to 

p(t) = ~[HcQD,p(t)]+V[T + Xm]p{t) (25a) 

= -tIKcqd - ih{T*X - W)m/2, p(t)} + V[Xn x + T + T]p{t), (25b) 
n 

= £p(t), (25c) 



where T> is defined in Eq. (|io|). Here T is an arbitrary complex numbeic3'c2l, while we are using T and X to represent 
respectively the quantities 7+ and X + in Eq. (^]) evaluated at zero temperature. 

Requiring that the ensemble average of the conditioned density matrix E[p c (t + dt)] — p(t + dt) satisfies the 
unconditional master equation ( p5|) leads to 

p 0c (t + dt) + p lc (t + dt) = (1 + dtC) Pc (t). (26) 

Here we have explicitly used the stochastic Ito calculusEl'H for the definition of time derivatives as p(t) — 
lim.dt^o[p{t + dt) — p(t)]/dt. This is in contrast to thq-dpfinition, p(t) = limdt— >n \p(t + dt/2) — pit — dt/2)]/dt, used in 



another stochastic calculus, the Stratonovich calculusE3S. Recall that Eq. ( |22|) indicates that E[dN c (t)]/dt equals to 
the average electron tunneling rate through the PC barrier. From Eq. (|[), the electron tunneling rates are D = \T\ 2 
when rii — and D' = \T + X\ 2 when n\ = 1. From Eq. (|22|) we thus have the correspondence 



Ti[Mi(dt)p c {t)Ml(d£)] = Tr{p c (i)[T* + n lX *][T + n lX ]}dt. (27) 



G 



Also, for Eq. (|26j) to reproduce the master equation ( p5b| ) we must have 

Mi(dt) = Vdl(Xni + T + T) (28) 
for some arbitrary complex number T . By inspection of Eq. ( p7|) we must have T = 0, so that 

p lc (i + ^) = J[T + X ni ]p c (t)dt. (29) 

Substituting Eq. (§|) into @ yields 

E[dN c (t)] = Tr[/5 lc (t + dt)} = [D + (£>' - D)(ni) c (i)]di, (30) 

where (ni) c (t) = Tr[ni/9 c (i)]. The remaining part, except the jump of Eq. (p9|), on the right hand side of Eq. ( p6| ) in 
time eft, corresponds to the effect of a measurement giving a null result on pjt): 

p 0c (t + dt) = p e (t)-dtLA[T + Xni]p e (t)-l[HcQD,Pc(t)]}, (31) 

where A is defined in Eq. (|l2|). The corresponding measurement operator is 

M Q {dt) = 1 - dt[(i/h)H C QD + (1/2)(T* + **m)(T + *m)]. (32) 

Finally, substituting Eqs. (|29]), (|3l|), ( pi| ) and ( |30| ) into Eq. (^), expanding and keeping the terms of first order in 
dt, we obtain the stochastic master equation, conditioned on the observed event in time dt: 



dp c {t) = dN c (t) 



J[T + X ni ] 

Plc(t) 



- 1 



where 



Pc {t) + dt i —A[T + Xm}p c (t) + V lc {t) Pc (t) + -[H C QD,Pc(t]] 



P lc (t) = D + (D'-D){m) c (t). 



(33) 



(34) 



Note that dN c (t), from Eq. (|30|), is of order dt. Hence terms proportional to dN c (t)dt are ignored in Eq. (|33j). Again 
averaging this equation over the observed stochastic pro cess by setting E[dN c (t)\ equal to its expected value, Eq. pCj), 
gives the unconditional, deterministic master equation ( 25a ). Eq. ( |33| ) is one of the main results in this paper. 

So far we have assumed perfect detection or efficient measurement. In this case, the stochastic master equation for 
the conditioned density matrix operator (|3^) is equivalent to the following stochastic Schodinger equation (SSE) for 
the conditioned state vector: 



d\Mt)) = 



dN c (t) 



T + Xm 



- 1 



- dt ( -TCcqd 



{T* +X*n 1 )(T+Xn 1 ) V lc (t) 



This equivalence can be easily verified using the stochastic Ito calcului 



d Pc (t) = d(\Mt))(Mt)\) 

= (d\Mt)))(Mt)\ + \Mt))d(Mt)\ 



(d\Mt)))(d(Mt)\), 



(35) 



(36) 



and keeping terms up to order dt. Since the evolution of the system can be described by a ket state vector, it is 
obvious that an efficient measurement or perfect detection preserves state purity if the initial state is a pure state. 
In this description of the SSE, the quantum average is now defined, for example, as (ni) c (t) — {ip c {t)\rii \ip c (t)). The 
unconditional density matrix operator is equivalent to the ensemble average of quantum trajectories generated by the 
SSE, p(t) — E[\^p c (t)\M c (t)\], provided that the initial density operator can be written as p(0) = IV-'c(O)) (V'c(0)|. 

The interpretation^ for the measured system state conditioned on the measurement, in terms of gain and loss of 
information, can be summarized and understood as follows. In order for the system to be continuously described 
by a state vector (rather than a general density matrix), it is necessary (and sufficient) to have maximal knowledge 
of its change of state. This requires perfect detection or efficient measurement, which recovers and contains all 
the information lost from the system to the reservoirs. If the detection is not perfect and some information about 
the system is untraceable, the evolution of the system can no longer be described by a pure state vector. For the 
extreme case of zero efficiency detection, the information (measurement results at the detector) carried away from the 
system to the reservoirs is (are) completely ignored, so that the stochastic master equation ( |33| ) after b eing averaged 
over all possible measurement records reduces to the unconditional, deterministic master equation (25a), leading 
to decoherence for the system. This interpretation highlights the fact that a density matrix operator description 
of a quantum state is only necessary when information is lost irretrievably. The purity-preserving, conditional State 
evolution for a pure initial state, and gradual purification for a non-pure initial state have been discussed in RefsJ13~E3't^l 
in the quantum diffusive limit. 
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IV. QUANTUM-DIFFUSIVE, CONDITIONAL MASTER EQUATION 







In this section, we extend the results obtained in previous section and derive the conditional master equation when 
the average electron tunneling current is very large compared to the extra change of the tunneling current due to the 
presence of the electron-in-thc dot closer to the PC. This limit is studied and called a "weakly coupling or responding 
detector" limit in Refs.BES. Here, on the other hand, we will refer to this case as quantum diffusion in contrast to 
the case of quantum jumps. In the quantum-diffusive limit, many electrons, (N > \(D' + D)/(D' — D)] 2 3> 1), pass 
through the PC before one can distinguish which dot is occupied. In addition, individual electrons tunneling through 
the PC are ignored and time averaging of the currents is performed. This allows electron counts, or accumulate 
electron number, to be considered as a continuous variable satisfying a Gaussian white noise distribution. In RefsJ13 : 
a set of Langevin equations for the random evolution of the CQD system density matrix elements conditioned on the 
detector output was presented, based only on basic physical reasoning. In this section, we show, exnlicitly, under the 
quantum-diffusive limit, that our microscopic approach reproduces^ 3 ] the rate gauations in Refs.ESo. 

In quantum optics, a measurement scheme known as homodyne detectionEilc3e3 is closely related to the measure- 
ment of the CQD system by a weakly responding PC detector. In both cases, there is a large parameter to allow 
the photocurrent or electron current to be approximated by a continuous function of time. We will follow closely the 
derivation of a smooth master equation for homodyne detection given in Ref.ca (sketched first by CarmichaeB) for 
the CQD system. 

There are two ideal parameters T and X for the CQD system. In the quantinHjjdjffusive limit, we assume |T| 3> \X\ 
which is consistent with the assumption, (D + D') {D 1 — D), made in RefsJl3 : E] for the weakly coupling or weakly 
responding PC detector. Consider the evolution of the system over the short-time interval [t, t+St). We relate the three 
parameters, X, T and St, in our problem as \X\ 2 St ~ e 3 / 2 , where e = (|<Y|/|T|) <C 1. This scaling is chosen so that in 
time St, the number of detections (electron counts) with dot 1 being unoccupied scales as SN ~ |T| 2 5t ~ (T 1 ! 2 ^> 1. 
However, the extra change in electron number detections due to the presence of the electron in dot 1 scales as 
\X\ 2 St ~ e 3 / 2 < 1. To be more specific, the average number of detections, following Eq. (^0|), up to order of e 1 / 2 is 



E[SN(t)} = \T\ 2 6t[l + 2e cos0 {ni) c {t)\, (37) 

where 9 is the relative phase between X and T. The variance in SN will be dominated by the Poisson statistics of the 
current eD = e\T\ 2 in time St. Since-the number of counts in time St is very large, the statistics will be approximately 
Gaussian. Indeed, it has been showncH that the statistics of SN are consistent with that of a Gaussian random variable 
of mean given by Eq. ( j37|) and the variance up to order of eT l l 2 is o 2 N — \T\ 2 St. The fluctuations a 2 N is necessarily as 
large as expressed here in order for the statistics of SN to be consistent with Gaussian statistics. Thus, SN can be 
approximately written as a continuous Gaussian random variableESw 

SN{t) = {\T\ 2 [l + 2e cos0(ni) c (t)] + \T\£(t)}St, (38) 

where £(i) is a Gaussian white noise characterized by 

E[£(t)]=0, E[£(t)i;(t')]=5(t-t'). (39) 



Here E denotes an ensemble average and S(t — t') is a delta function. In stochastic calculusLnE 3 ], £(t)dt = dW[t) is 
known as the infinitesimal Wiener increment. In Eq. (p8|), the accuracy in each term is only as great as the highest 
order expression in e 1 / 2 . But it is sufficient for the discussions below. 

Although the conditional master equation ( |33| ) requires that dN c (t) to be a point process, it is possible, in the 
quantum-diffusive limit, to simply replace dN c (t) by the continuous random variable SN c (t), Eq. (|38|). This is because 
each jump is infinitesimal, so the effect of many jumps is approximately equal to the effect of one jump scaled by 
the number of jumps. This can be justified more rigorously as in Ref.Ell. Finally, expanding Eq. ( p3| ) in power of 
e, substituting dN c (t) — > SN c (t), keeping only the terms up to the order e 3 / 2 , and letting St — > dt, we obtain the 
conditional master equation 

f>S) = -^[H C QD,Pc{t)}+V[T + Xm]p c (t) 

[T*Xn lPc (t)+X*T Pc (t)ni -2Re(T*X){n 1 } c (t)p c (t)\ . (40) 

Thus the quantum-jump evolution of Eq. (33) has been replaced by quantum-diffusive evolution, Eq. (^0[). Following 
the same reasoning in obtaining the SSE (35) for the case of quantum-jump process, we find the quantum-diffusive, 
conditional master equation (EQ) is equivalent to the following diffusive SSE: 
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d\Mt)) 



dt 



-7~Lcqd 



\X\ 



[m - 2ni(m) c (t) + (ni) 2 c (t)} -iIm(T*X)ni 



+ f(t)*-=T {T*Xm - X*T (m) c (t)} 



(41) 



that 



This equivalence can be verified using Eq. (36) and keeping terms up to oder dt. Note however in this case 
terms of order £(t)dt are to be regarded as the same order as dt, but [£,(t)dt] 2 = [dW(t)] 2 — dt. 

Our conditional master equation by its derivation is formulated. in, terms of Ito calculus, while the stochastic rate 
equations in RefsJl3'E3 are written in a Stratonovich calculus formEjO. In contrast to the Stratonovich form of the 
rate equations, it is easy to see that the ensemble average evolution of our conditional master equation ( [40| ) reproduces 
the unconditional master equation (25a) by simply eliminating the white noise term using Eq. (|39|). To show that 
our quantum-diffusive, conditional stochas tic m aster equation (^p| ) r e produces the non-linear Langevin rate equations 
obtained semi-phenomenologically in RefsJl3E3, we evaluate Eq. (ffO[) in the same basis as for Eq. (13) and obtain: 



f>aa(t) = ifl[pab(t) - pba(t)] ~ 
Pab(t) = i£Pab(t) + itt[Paa{t) 



/$r d Paa(t)p bb (t)Z{t) 

Pbb(t)} - T d p ab {t) + 



/ ^'dPab{t)[p aa {t)- f>bb(t)]Z(t) 



(42a) 
(42b) 



In obtaining Eq. (|42j), we have made the assumption of real tunneling amplitudes as in Refs.li3ii§E!l in order to be able 
to compare the results directly. We have also set X — —y/2Td- Again, the ensemble average of Eq. (42) by eliminating 
the white noise tetjas. reduces to Eq. (p"3|). To further demonstrate the equivalence, we translate the stochastic rate 
equations of Refs.E3'c£l into Ito formalism and compare them to .. E |q. (|42|). This is carried out in Appendix |a|. Indeed, 
Eq. (E3) is equivalent to the Langevin rate equations in Refs.li3iH3 for the "ideal detector" . 



V. ANALYTICAL RESULTS FOR CONDITIONAL DYNAMICS 



To analyze the dynamics of a two-state system, such as the CQD system considered here, one can represent the 
system density matrix elements in terms of Bloch sphere variables. The Bloch sphere representation is equivalent to 
that of the rate equations. However, some physical insights into the dynamics of the system can sometimes be more 
easily visualized in this representation. Denoting the averages of the operators <r x , cr y , a z by x, y, z respectively, the 
density matrix operator for the CQD system can be expressed in terms of the Bloch sphere vector (x,y, z) as: 



p(t) = [I + x(t)cr x + y(t)a y + z{t)a z ]/2 

If 1 + z(t) x(t) - iy{t) 
~ 2 \x(t)+iy(t) l-z(t) 

where the operator /, at, are defined using the fermion operators for the two dots: 

I = c|c2 + c\c\, 

(Jy = — «4 C 1 + « C 1 C 2> 

0- 2 = 4 C 2 - C 1 C 1- 



(43a) 
(43b) 



(44a) 
(44b) 
(44c) 
(44d) 



It is easy to see that Trp(i) = 1, J is a unit operator, and <Ji defined above satisfies the properties of Pauli matrices. In 
this representation, the variable z(t) represents the population difference between the two dots. Especially, z(t) = 1 
and z(t) = —1 indicate that the electron is localized in dot 2 and dot 1 respectively. The value z(t) = corresponds 
to an equal probability fo r the electron to be in each dot. 

The master equations (25a), ( |40| ) and(|33|), can be written as a set of coupled stochastic differential equations in 
terms of the Bloch s pher e variable s. Fo r simplicity, in this section we assume that the tunneling amplitudes are real. 
By substituting Eq. (43a) into Eq. (25a), and collecting and equating the coefficients in front of a x , a y , a z respectively, 



the unconditional master equation under the assumption of real tunneling amplitudes is equivalent to the following 
equations: 



dt \ y(t) ) \ e -T d 
dz(t) 



x(t) 

y(t) 



o 

-2flz(t) 



dt 



= 2n y (t). 



(45a) 
(45b) 
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Similarly for the quantum-diffusive, conditional master equation JIc|), we obtain 



dx c (t) 

dt 
dy c {t) 

dt 
dz c (t) 

dt 



~ey c {t) - T d x c (t) - ^2T^z c (t)x c (t)^(t), 

ex c (t) - 2il z c (t) - T d y c (t) - ^2T~ d z c (t)y c (t)£{t), 

2Qy c (t) + y/2~F d [l~z 2 c (t)] £(t). 



(46a) 
(46b) 
(46c) 



Again the c-subscript is to emphasize that these variables refer to the conditional state. It is trivial to see that Eq. 
( p6| ) averaged over the white noise reduces to Eq. fl45|), provided that E[x c (t)] = x(t) as well as similar replacements 
are performed for y c (t) and z c {t). The analogous calculation can be carried out for the quantum-jump, conditional 
master equation (j33l). We obtain 



(47a) 
(47b) 
(47c) 



'/•'■•■(' ) = ( -eyS) (D \ D) z c (t)x c (t)) ~ dN c (t) ( x&) 2Td (D ' D)Zc[t) 



dy c (t) = dt ( ex c {t) - 2flz c {t) - ^ D \ c {t)y c {t) ) - dN c (t) ( y c (t) 



2D + (D' -D)[l-z c (t)}, 

2T d - (£>' - D)z c {t) 



dz c (t) = dt ( 2Qy c (t)dt + ^— ^ [l - z 2 c {t)] ) - dN c (t) 



2D +{D'-D) \l-z c {t)] 
(D>-D){l~z 2 c (t)} 



2D +(D'-D) [l-z c (t)\ 



As expected, by using Eq. (|30D, the ensemble average of Eq. (|47]) also reduces to the unconditional equation (|45|). 

Next we calculate the localization rate, at which the electron becomes localized in one of the two dots due to the 
measurement, using Eqs. @) and @). Obviously, the stochastic, conditional differential equations provide more 
information than the unconditional ones do. In the unconditional case Eq. (f45|), the average population difference 
z(t) between the dots is a constant of motion ([dz(t)/dt] = 0) when = 0. However, if the present model indeed 
describes a measurement of n% = c\c\ (in other words the position of the electron in the dots), then in the absence 
of tunneling Q = 0, we would expect to see the conditional state become localized in one of the two dots, i.e., either 
2=1 or z = — 1. Indeed, for = 0, we can see from the conditional equations (46c) and (47c) that z c (t) — ±1 are 
fixed points. We can take into account both fixed points by considering z 2 c (t). Hence it is sensible to take the ensemble 



average E[z^(t)] and find the rate at which this deterministic quantity approaches one. applying ltd calculusE3H to 
the stochastic variable z„(t ), we have d(z^) — 2z c dz c + dz c dz c . Let us first consider the case for the quantum-jump 
equations. Using Eqs. (47c) and j30| ) and the fact that dN^{t) = dN c (t), we find that 



E[dz 2 c {t)] = E 



{D'-D?[l-zm 2 
4D + 2(D'-D)[l-z c {t)} 



dt. 



(48) 



If the system starts in a state which has an equal probability for the electron to be in each dot then z c (0) = z(0) = 0. 
In this case, the ensemble average variable z(t) would remain to be zero since [dz(t)/dt] = when f2 = 0. However 
if we average z 2 over many quantum trajectories with this initial condition then we find from Eq. ( f48| ) that for short 
times (by setting z c (t),z 2 (t) m 0) 



E[z 2 (5t)} 



(D' - D) 2 
2(D' + D) 



That is to say, the system tends toward a definite state (with z c = ±1 so = 1) at an initial rate of 



1 ump 
Tloc = 



(D' — D) 2 _ (V^ + VD) 
2(D' + D) ~ (£>' + D) 



(49) 



(50) 



Similarly for the case of quantum diffusion, using Eqs. (kl7d) and (B9J) and the fact that [£(t)dt} 2 = [dW{t)} 2 = dt, 



we find E[dz 2 (t)] = E [21^ [1 — z^(i)] 2 ] dt. Applying the same reasoning for obtaining Eq. ([49]), we find E[z 2 (St)] 



2T d St = 7^'* St. This implies that the localization p.te in this case is 7^ = 2T d . This is consistent with the result 
of localization time, t\ oc ~ (l/7io?)j f° un d m Ref.ta in the quantum-diffusive case. As expected, Eq. (50) in the 
quantum-diffusive limit, T 3> X or (D + D') ^> (D' — D), reduces to 7{"™ p — > 2T d = 7^. The rate of localization is 
a direct indication of the quality of measurement. It is necessarily as large as the decoherence rate since a successful 
measurement distinguishing the location of the electron on the two dots would destroy any coherence between them. 
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The above localization rates are related to the signal-to-noise ratio for the measurement and can be obtained 
intuitively as follows. Consider the electron with equal likelihood in either dot so that z c (0) = z(0) = 0. For the case 
of quantum diffusion, the electron tunneling current through the PC obeys Gaussian statistics. Recall in Sec. jrV] that 
the mean of the probability distribution of the number of electron detections through the PC is given by Eq. (|37| ) and 
its variance takes the form af^ — T 2 5t in time St. If the electron is in dot 1, then the rate of electrons passing through 
the PC is T 2 + 2TX; if it is in dot 2, then the rate is T 2 . One may define the width of the probability distribution 
as the distance from the mean when the distribution falls to e" 1 of its maximum value. For a Gaussian distribution, 
the square of the width is twice the variance. The above two probability distributions will begin to be distinguishable 
when the difference in the means of the two distributions is of order the square root of the sum of twice the variances 
(square of the widths) at time r. That is, when 



[(T 2 + 2TX)t - T 2 t] ~ y/2T 2 r + 2T 2 r. (51) 

Solving this for r gives a characteristic rate: r _1 ~ X 2 = 2T,i- This is just the 7^ discussed above. For the case of 
quantum jumps, the statistics of the electron counts through the PC can be approximated by Poisson statistics. For 
a Poisson process at rate 1Z, the probability for N events to occur in time t is 

p(tf;t) = £*£e-« (52) 
The mean and variance of this distribution are equal and given by E[N] — Vslt(N) = IZt. In the quantum-jump case 



from Eq. (30), if the electron is in dot 1 then the rate of electrons passing through the PC is D' . If the electron is in 
dot 2, then the rate is just D. Requiring the difference in means of the two probability distributions, p(N,r), being 
of order the square root of the sum of twice the variances at time r yields: 



(D't - Dt) ~ V2D't + 2Dt. (53) 

Solving this for t _1 yields a characteristic^iate^hich is the same as 7^™ p defined in Ea^IEOL 

A similar conclusion is reached in Refs.EHO'El. The measurement time, i ms , in RefsE3'c3E3 is roughly the inverse 
of the localization rate given here. However, there the condition for being able to distinguishJiie two probability 
distribution is slightly different from the condition discussed here. The measurement timeEZrEilO is denoted as the 
time at which the separation in the means of the two distributions is larger than the sum of the widths, i.e., the sum 
of the square roots of twice the individual variance rather than the square root of the sum of twice the variance. If 



this condition is applied here, instead of Eqs. (|5l|) and (53), we have 



[(T 2 + 2TX)t ms - TH ma ] > ^/2TH ms + ^2TH ms (54) 
for the quantum-diffusive case, and 



(£>' t ms - D t ms ) > ^2D't ms + V2Dt ms (55) 

for the quantum-jump case. We find from Eqs. (|54|) and ( |55| ) that the inverse of the measurement time f ms is the 
same for both quantum-diffusive and quantum-jump cases, and is equal to the decoherence rate: 

x 2 (Vd>-Vd) 2 _ t 

*ms = — = 2 = T d = T d , (56) 

where = (l/T^) is the decoherence time. This is in agreement with the result in Refs. Our condition shows, on 
the other hand, the different localization rates for the quantum-jump and quantum-diffusive cases. This is consistent 
with the initial rates obtained from the ensemble average of Bloch variable, E\^W^L 

There is another time scale denoted as mixing time, i m i x , discussed in Refs.EZrEj'L 2 ]. It is the time after which the 
information about the initial quantum state of the CQDs is lost due to the measurement-induced transition. This 
transition arises because of the non-zero coupling Q term in the CQD Hamiltonian, which does not commute with 
the occupation number operator of dot 1 (the measured quantity) and thus mixes the two possible states of the CQD 
system. Below we estimate the mixing time using the differential equations for the Bloch variables. It is expected 
that effective and successful quantum measurements require t m i X 3> t ms ~ t\ oc ~ r^. In other words, the readout 
should be achieved long before the information about the measured initial quantum state is lost. In terms of different 
characteristic rates, we have, in this case, the relation: Td ~ 710c ~ £ms ^ 7mix, where 7 m i X = (l/imix) represents the 
mixing rate. For finite O, the rate at which the variables x(t) and y(t) relax can be found from the real part of the 
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eigenvalues of the matrix in the first term on the right hand side of Eq. ( J45a| ) . This gives the decay (decoherence) rate 
Trf for the off-diagonal variables, x(t) and y(t). The variable z{t) = represents an equal probability for the electron 
in the CQDs to be in each dot. Hence the rate at which the variable z(t) relaxes to zero corresponds to the mixing 
rateE3, 7 m i X . Under the assumption of Yd 3> 7mix for effective measurements, the variables x(t) and y(t) therefore 
relax at a rate much faster tha n that of the v ariab le z(t). As a result, it is valid to substitute the steady-state value 
of y(t) obtained from Eq. (45a) into z{t) Eq. (45t) to find the mixing rate. Consequently, we obtain 



dz(t) _ 4n 2 r d 

~dT ~ ~T 2 d + e 2 



z(t) = -7mixZ(t). 



(57) 



It is easy to sec that the mixing rate Eq. ( 
assumption ^> 7 m ; x yields, from Eq_. ( 



the condition: ft <C ( \/T d + e. 
agreement with the result found in Ref.Ea under the similar required condition! 



vanishes as SI — > 0. Finally, the self-consistent requirement for the 

2). 



The mixing rate Eq. (B7J) is in 



VI. CONCLUSION 



We have obtained the unconditional master equation for the CQD system, taking into account the effect of finite- 
temperature of the PC reservoirs under the weak system-environment coupling and Markovian approximations. We 
have also presented a quantum trajectory approach to derive, for both quantum-jump and quantum-diffusive cases, 
the zero-temperature conditional master equations. These conditional master equations describe the evolution of the 
measured CQD system, conditioned on a particular realization of the measured current. We have found in both cases 
that the dynamics of the CQD system can be described by the SSEs for its conditional state vector provided that the 
information carried away from the system by the PC reservoirs can be recovered by perfect measurement detection. 
Furthermore, we have analyzed for both cases the localization rates at which the electron becomes localized in one 
of the two dots when $1 = 0. We hai[e_&h.own that the localization time discussed here is slightly different from the 
measurement time defined in Refs.EDOH3. The mixing rate at which the two possible states-pf the CQDs become 
mixed when f2 ^ has been calculated as well and found in agreement with ihe result in Ref.E2l. 

We thank M. Biittiker and A. M. Martin for bringing our attention to RefEj. HSG is grateful to useful discussions 
with A. N. Korotkov, D. V. Averin, K. K. Likharev, D. P. DiVincenzo, J. R. Friedman and V. Sverdlov. 



APPENDIX A: EQUIVALENCE OF STOCHASTIC RATE EQUATIONS IN DIFFERENT CALCULUS 

In this Appendix, we translate the stochastic rate equations of Refs. EE , written in terms of Stratonovich calculus, 
into Ito calculus formalismHEj. Although the translation was sketched and the result was stated in Ref.u, for 
completeness, we fill in the calculation steps using our notation here. Eq. (11) and (12) of Ref.liS in Stratonovich 
calculus formalism are rewritten in terms of our notation as follows: 



p bb (t) = *n[p6a(t) - Pab{t)] + 4^ p bb {t)p a<1 {t) y-VSiT d [p aa {t) - p bb {t)] + \J yf (*) J , (Al) 

Pba(t) = ~isp ba (t) + iQ[p bb (t) - paa{t)] 

r ^ „m ( /TTT^r ^ ^ , [ft 



+ 2^/ ^ [ P aa(t) ~ Pbb (t)} i-VS^iPaait) - Pbb(t)} + yj y^)J Pba(t), (A2) 

where we have substituted the notation used in R(i£t3 to H/fi — > fl, e/fi — > e, and expressions for 1Z and AI in 
terms of Td and Sj using Eqs. (10) and (2) of RefJla. Specifically, we have set AI = —2y/SiTd- In addition, the 
white noise £(t) in RefJla has spectral density S% = Si, which implies E[£(t)£(t')] = (Si/2)S(t — t'), different from our 



definition, Eq. (|39|). Hence, the replacement £(fL-j-> y/ Si/2 £(t) has been employed. Moreover, since an ideal detector 
is assumed, 7^ is set to zero for Eq. (12) of Ref.ES. Note finally that the electron operator indices in the CQD should 
be interchanged. Foe-,example, the electron annihilation operator C2 in Ref.EEl should be c\ in our notations. As a 
result, pn{t) in Ref.E3 is rewritten as p bb (t), and pig(t) as p ba (t) here. 



As pointed out in RefB, to translate Eqs. ( |Al| ) and ( |A2| ) into Ito formalism, one needs to add the terrrB 
(F/2)(dF/dpij) for each rate equation f>ij(t), where F is the factor before £(t) in each equation respectively. Note 
that the factor Si /2 appearing in front of the term needed to be added for the translation in Ref .EJ is set to 1 here. 
This is because of the different definitions of the stochastic white noise variables £(t) in both cases, discussed above. 
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To be more specific, F — y/8Td pbb{t)p aa (t) for Eq. (Al) in our notations. By using the relation p aa (t) = 1 — Pbb(t), 
it is easy to find the derivative dF/dpbb = V^Td [p a a(t) — pbb(t)]- As a consequence, the term needed to be added to 
Eqs. (CO) is 



4FdPbb(t)p aa (t)[paa(t) - Pbb(t)}, 



(A3) 



which exactly ca ncel s the first term inside t he bi g parenthesis in the second term of Eq. (Al). Hence the resultant 
equation for Eq. (Al) in ltd form is just Eq. ( |42a| ) with an overall minus sign in front of it (pbb(t) — —f>aa(t))- As for 
Eq. ( |A2| ), it is easy to find that F = \/2T c i [p aa {t)^= Pbb(t)]Pba(t)- In order to carry out the derivative with respect to 
Pba{t)7o^e needs the expression of Eq. (8) of Ref.Ea to relate diagonal elements, pbb(t) and p aa (t), to pb a {t)- We then 
obtain 



dF(t) 



2[p aa (t) - pbbjt)} 2 - 1 

Paa(t) - Pbb{t) 



dpba(t) 

Thus the terms needed to be added into Eq. ( |A2| ) are 

2T d [p aa (t) ~ Pbb(t)} 2 Pba{t) ~ T d p ba {t). 



(A4) 



(A5) 



The first term in Eq. (A5) exactly cancels the term with the square bracket inside the big parenthesis in the last term 
of Eq. (|A2|). Therefore the resultant equation for Eq. ( A2) in ltd form is equal to the complex conjugate of Eq. (12b) 
(pba(t) = P* a b(t))- This completes our demonstration of the equivalence. 
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